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We examine flat models containing a dark matter component and an arbitrary dark energy com- 
ponent, subject only to the constraint that the dark energy satisfies the weak energy condition. We 
determine the constraints that these conditions place on the evolution of the Hubble parameter with 
redshift, H(z), and on the scaling of the coordinate distance with redshift, r(z). Observational con- 
straints on H(z) are used to derive an upper bound on the current matter density. We demonstrate 
how the weak energy condition constrains fitting functions for r(z). 
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I. INTRODUCTION 

Observational evidence [3, Q indicates that roughly 
70% of the energy density in the universe is in the form 
of an exotic, negative-pressure component, dubbed dark 
energy. (See Ref. fH for a recent review). If pde and 
Pde are the density and pressure, respectively, of the 
dark energy, then the dark energy can be characterized 
by the equation of state parameter w, defined by 



w 



Pde/p 



DE- 



(1) 



Although the simplest possibility for the dark energy is 
a cosmological constant, which has w = — 1, many other 
possibilities have been proposed, including an evolving 
scalar field (quintessence) [3, H, H 0, H], a scalar field 
with a non-standard kinetic term (fc-essence) 0, [13, HH, 

m, m, q m m, El , or sim P i y an arbitrary barotropic 

fluid with a pre-determined form for p(p), such as the 
Chaplygin gas and its various generalizations [l|| I20I 

Lacking a definite model for the dark energy (aside 
from the perennial favorite cosmological constant), it is 
interesting to determine what can be derived from fairly 
general assumptions about the nature of the dark energy. 
In particular, a plausible assumption about the dark en- 
ergy is that it obeys the weak energy condition (WEC). 
If T^ v is the energy- momentum tensor of the dark energy, 
then the weak energy condition states that 



TUff > 0, 



(2) 



where is any timelike vector. In a Friedman- 
Robertson- Walker universe, this reduces to a condition 
on the density and pressure: 



and 



P>0, 



P + P > o. 



(3) 



(4) 



Although models in which the dark energy violates the 
WEC are not inconsistent with current observations (as 
first pointed out by Caldwell [26|]), there are goo d reasons 
to believe that the WEC is satisfied [13, H, [H [H . 



The WEC has already been used previously to con- 
strain the expansion history of the universe [3l|, [H, [H, 
[H, [H, [H, [33| ■ These previous studies, however, all ap- 
plied the WEC constraints (equations [3] and [4j to the 
total cosmological fluid. Similarly, Schuecker et al. [HI 
applied several other energy conditions to the total cos- 
mological fluid. In this paper, instead, we assume a 
model consisting of matter plus a fluid obeying the WEC, 
and we determine the corresponding constraints that this 
places on the expansion history. 

In the next section, we derive the limits that can be 
placed on H(z) and r(z) from the WEC. In Sec. Ill, we 
apply these constraints to analyses of the observations. 
Our results are discussed in Sec. IV. 



II. CONSEQUENCES OF THE WEAK ENERGY 
CONDITION 



We assume a flat Friedman-Robertson- Walker model, 
containing a pressureless matter component and a dark 
energy component obeying the WEC. The radiation con- 
tribution to the density at late times is negligible and can 
be neglected. The matter density includes both baryonic 
and dark matter, and it scales with redshift z as 



Pm = Pmo(1 + zf 



(5) 



where the subscript will refer throughout to present- 
day values. The WEC imposes two conditions on the 
dark energy density: first, that 



Pde > 0, 



(6) 



at any redshift (a restatement of equation[3]) and second, 
that 



dppE 
dz 



>0, 



(7) 



which is a consequence of equation (j4]) . 

Now consider the consequences of these constraints for 
the redshift-dependent Hubble parameter H(z), defined 
by 



H{zf 



8vrG 



(Pm + Pde)- 



It is convenient to work in terms of the present-day value 
Hubble parameter, H , and to use the critical density pc 
defined by Hq = (8irG/3)pc- As usual, we take Qm = 
Pmo/pc and &de = Pdeo/pc, and our assumption of a 
flat universe gives Qm + &>de = 1- Our fi's will always 
refer to present-day quantities; we omit the subscript 
in this case for simplicity. Then equation ([8]) can be 
rewritten as 

H{zf = n M (i + zf + (i - n M )( PDE (z)/ PDE0 ), (9) 

where we have defined H(z) = H(z)/H Q . The WEC 
forces (pde(z)/pdeo) > 1 for z > 0, so equation ([9|) 
becomes 



R{zf > fl M {l + z) 3 +1 - fi 



(10) 



Going back to equation ^ and taking the derivative with 
respect to z (which we denote throughout with a prime) 
gives: 

2H(z)H'(z) = 3fi M (l + z? + [(1 - ^mVpdeo^P-- 

az 



Now the WEC implies that dpoE/dz > 0, so we get 



(11) 




FIG. 1: Upper bound on the present-day value of Qm from 
the weak energy condition applied to the dark energy, using 
the H(z) values in Ref. for H = 73 ± 3 km s" 1 Mpc -1 . 



H{z)H\z) > -n M (l + z) 2 . 



(12) 



Equations (|T0[) and ([12]) give the constraints that the 
WEC for the dark energy places on the redshift- 
dependent Hubble parameter. Note that these con- 
straints are not independent. The limit in equation (fT2"]) . 
together with H(0) = 1, implies the limit given in equa- 
tion (fl~0|) . but the converse is not true. Although equation 
(fT2|) gives the stronger limit, we include both limits be- 
cause current data can provide some estimates for H(z) 
(as in the next section), but are too poor to provide any 
limits on H'{z). By construction, the standard ACDM 
model saturates both limits. Equation (|12p was previ- 
ously introduced by Sahni and Starobinsky in the con- 
text of quintessence models [33| and by Boisseau et al. 
in the examination of scalar-tensor models [40| . 

Now consider the coordinate distance r(z), defined by 



r{z) 



dz' 



(13) 



The coordinate distance is important because it is di- 
rectly related to the luminosity distance, d E , through 
d.L = c(l + z)r(z) / Hq, and it is d^ which is measured 
in supernova redshift surveys. Hence, a considerable ef- 
fort has been put into desig ning pa rametrizations for r(z ) 
to fit to the supernova data filt S lil lij. liol. Eol WA ltj| . 
The purpose of many of these investigations is to go from 
a best-fit form for r(z) to the potential for an underlying 
quintessence model. 

The derivative of equation (fTS)) gives 



r'(z) 



1 



H(zY 



(14) 



which, when combined with equation (jlOp . yields 

r'(z) < H^[n M (i + zf + (i - n M )}- 1/2 . 

Similarly, equation (fl"2)) yields 



^4>|^m(1 + .) 2 . 
r {zy 2 



(15) 



(16) 



Equations (|15|) and (fl6|) give the dark-energy WEC con- 
straint on the coordinate distance. As for the limits on 
H(z), equation (fl"6|) implies (jX5j) , but the converse is not 
true. 



III. COMPARISON WITH OBSERVATIONS 

Consider first our limits on the evolution of H(z). 
Estimates of H(z) were derived by Simon, Verde, and 
Jimenez (49j using passively evolving galaxies; these lim- 
its have been used to constrain cosmological parameters 
in dark energy models [5(| • A second approach to deriv- 
ing H(z) based on the Supernova data has been explored 
m Refs. [IH and We will use the results of [49| in 
our discussion here. 

Equation (|10[) can be rewritten as 



n M < 



H(zf - 1 



(1 



- 1 



(17) 



Thus, a single value of H (z) can provide an upper bound 
on f^M- Since our bound applies to any dark energy 
model satisfying the WEC, we do not attempt to fit any 
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FIG. 2: As Fig. 1, for H = 68 ± 4 km s" 1 Mpc -1 . 



particular model (as was done in Ref. |50|): rather, we 
calculate the upper bound individually for each H(z) 
measurement. The upper bound on fijy depends on 
H(z) — H(z)/H , so our results will naturally be sen- 
sitive to the value of Ho- Following Ref. |50| we consider 
two priors for Hq: Hq = 73 ± 3 km s _1 Mpc -1 from 
WMAP [H, and H = 68 ± 4 km s" 1 Mpc" 1 from the 
median statistics analysis in Ref. (54J. Our results are 
shown in Figs. 1 and 2 (where all error bars are 1-sigma). 
For z = 1.53, where the error bars are smallest, we get 
the tightest upper bound on Qm- < 0.18 ± 0.05 for 
H = 73 ± 3 km s" 1 Mpc" 1 , and Q M < 0.21 ± 0.06 
for Ho = 68 ± 4 km s _1 Mpc -1 . (For a very different 
approach, see Ref. [55j|). 

Next we consider the consequences of the limits on 
the coordinate distance r(z) given by equations (fl~5f and 
(fT!))) . The quantity observers actually measure is the ap- 
parent magnitude m(z), given by 

m(z) = M + 5Log w {D L (z)), (18) 

where Dl is the Hubble- free luminosity distance, 

D L (z) = (l + z)J o dz'j^, (19) 

and M. is the magnitude zero-point offset, which depends 
on the absolute magnitude M as 

M = M + 5Log 1Q ( \ + 25 = M - 5Log w h + 42.38. 
\ Mpc J 

(20) 

The distance modulus given by the SNIa data is defined 
as 



FIG. 3: The region in the wi, wi plane allowed by the weak 
energy condition for the parametrization of r(z) given in 
equation (|24[1 . The region between the two sets of curves 
is allowed. Right-hand boundaries are given for £7 a/ = 0.2 
(dashed), Qm — 0.3 (dotted) and Qm = 0.4 (dot-dash). Left- 
hand boundary (solid) is independent of Qm ■ The two ellipses 
are the la and 2a contours obtained by fitting equation (|22|) 
to the supernova data without the weak energy constraint. 



where fio = 42.38 — bLogioh. 

As an example, we now choose a representative fitting 
function, translate it into a fitting function for r(z), and 
then determine how the WEC constrains the parame- 
ters of the fitting function. In this paper, we use the 
fitting function for /x, first introduced by Padmanabhan 
and Choudhury (47|, given by 



M/it = Mo + 5Log w 



z(l + W\Z) 



(l + w 2 z) 



(22) 



where /io? w i an d w 2 are the three independent fitting 
parameters. By comparing equations (21) and (22), we 
can write 



D L = 



so that 



r(z) 



1 



z(l + W\Z) 
(1 + W2Z) 



z(l + w\z) 



Ho (I + z){l + w 2 z) 



(23) 



(24) 



fi(z) = m(z) - M = 5Log w (D L (z)) + /i 0) 



(21) 



Since we have an analytic form for r(z), we use equation 
(fl"6|) alone to determine the values of w\ and W2 which vi- 
olate the WEC; any r(z) which satisfies equation (TIT)]) will 
automatically satisfy equation (|15p . We require equation 
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([T6|) to be satisfied for z < 2, the range over which the 
supernova data extend. 

The allowed region for w\ and u>2 is displayed in Fig. 
3. As expected, the constraints are tighter for larger val- 
ues of Omj but rather surprisingly, the excluded region 
in parameter space is rather insensitive to the assumed 
value of (and a significant region of parameter space 
is excluded even in the limit SIm — > 0). Also in Fig. 3, we 
display the confidence regions for w\ and W2 from the su- 
pernova data, with no assumptions about the equation of 
state for the dark energy. We use 60 Essence supernovae 
[56[, 57 SNLS supernovae [13] and 45 nearby supernovae. 
We have also included the new data release of 30 SNe la 
detected by HST and classified as the Gold sample by 
Riess et al. [Hj]. The combined data set can be found 
in Ref. [5{|. Clearly, the best-fit values for w\ and W2 
lie slightly outside of the allowed region for Qjm > 0.3. 
Not too much should be read into this: given that ACDM 
models provide a good fit to all current observations, and 
such models saturate our bounds, we would expect the 
best fit parameters to lie near the boundary of the ex- 
cluded region. The important point is that Fig. 3 shows 
how one can use our constraints on r(z) to eliminate re- 
gions of parameter space for which the dark energy vio- 
lates the WEC. 

Of course, it is also possible to parametrize w as a 
function of z and use this parametrization to derive a 
form for r{z) (see, e.g., Ref. (HJ); m this case the WEC 
is trivially satisfied as long as the parametrization for 
w(z) forces 1 + w(z) to be nonnegative. 



combined with observations of H (z) to put upper bounds 
on Qm- While the scatter in these estimated upper lim- 
its is large, as are the errors, the important point is that 
these are generic upper limits, independent of the na- 
ture of the dark energy (as long as it satisfies the weak 
energy condition). Improved measurements, particularly 
of H(z), will strongly improve this upper bound. The 
constraints on r(z) do not provide similarly useful lim- 
its, but they can be applied to any parametrization of 
r(z) to eliminate in advance any regions of parameter 
space in which the dark energy violates the weak energy 
condition. 

What happens if the WEC is violated by the dark 
energy? If one allows for arbitrary evolution, then 
there are clearly no constraints on r(z) and H(z). An 
intermediate case, which provides weaker limits than 
the ones we have discussed, is when the WEC applies 
to the total fluid (matter and dark energy together) 
[3l|,[33,[3i,[33,[3i,[3a[33]- In this case, for example, there 
is no limit corresponding to equation (flT))) . while the limit 
corresponding to equation (|12p becomes H{z)H' (z) > 0. 
Thus, this version of the WEC provides no bound on FIm, 
although it does constrain the evolution of H(z). Apply- 
ing the WEC to the total fluid also yields constraints on 
r(z); these are discussed in detail in Ref. [Hj]. These 
constraints are weaker but more general than the con- 
straints we have obtained by applying the WEC to the 
dark energy alone. 



IV. DISCUSSION 
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